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Abstract 

This paper develops an analytical relation on probability to ob- 
tain a condition on stability of steady states of open systems which 
represents the condition of quasi-ergodicity for the open systems too. 
It is proved that the result obtained is consequence of the maximum 
entropy generation principle. 
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1 Introduction 

In phenomena out of equilibrium irreversibility manifests itself because the 
fluctuations of the physical quantities, which bring the system apparently out 
of stationarity, occur symmetrically about their average values (Gallavotti 
2006). The e-stady state definition allows us to obtain that for certain fluctu- 
ations the probability of occurrence follows a universal law and the frequency 
of occurrence is controlled by a quantity that has been related to the entropy 
generation. Moreover, this last quantity has a purely mechanical interpreta- 
tion which is related to the the ergodic hypothesis which proposed that an 
isolated system evolves in time visiting all possible microscopic states. More- 
over, considering that the open system is a system with perfect accessibility 
represented as a probability space in which is defined a PA-measure and a 
Borel function process the ergodic hypothesis itself is a consequence of the 
£-steady state definition because of the Hypothesis U\ introduced in Section 
4. 



The principle of maximum entropy generation represents the macroscopic 
effect of these theories, which, conversely, are its statistical interpretation. As 
a consequence of this principle it has been proved the irreversible ergodicity 
(Lucia 2008). As H. Primas has pointed out (Primas 1999): 'The impor- 
tance of the ergodic theorem lies in the fact that in most applications we 
see only a single individual trajectory that is a particular realization of the 
thermostatic process. The Kolmogorov's theory of stochastic processes refers 
to equivalence classes of functions and the Birkoff 's ergodic theorem provides 
a crucial link between the ensamble description and the individual descrip- 
tion of chaotic phenomena.'. The modern concept of subjective probabilities 
imphes a behaviour based on the Boolean logic (Primas 1999). The result 
obtained represents also the relation between the probability measurement in 
the probability space and real statistical facts. Following H. Primas (Primas 
1999) the set-theoretical representation of the Boolean algebra in terms of 
a Kolmogorov probability space is useful because it allows to relate a dy- 
namic in terms of a probabilistic density and all known context-independent 
physical laws are formulated in terms of pure states (Primas 1999). 

By entropy generation it has been introduced the global thermodynamic 
effect of the action of the forces finding a global expression which links the 
microscopical values to the macroscopical thermodynamic quantities (Lu- 
cia 2008). On the other hand, it has been argued that the information of 
a complex system may be incomplete because the dynamics of the system 
can be only partially known consequently to its not complete accessibility 
(Wang 2004). In addition the PA-measure in the probability space is not 
yet univocally defined, and it is difficulty to develop a statistical thermo- 
dynamics for an irreversible system. In 2009 the integral expression of the 
PA-measure in the probability space, useful in statistical thermodynamics 
of the complex and irreversible systems, has been obtained (Lucia 2009). To 
do it the basis of the incomplete information on complex system was devel- 
oped, underlying that in non-equilibrium transformation the volume of the 
phase space contracts; consequently, entropy generation were related to the 
incomplete information theory by using the Wang development of the Tsallis 
theory. This result underlines that the relevant physical quantities for the 
stochastic analysis of the irreversibility are the probability of the state in 
the phase space, the Hamiltonian valued at the end points of trajectory, and 
time (Lucia 2009). 

But, the Hamiltonian valued at the end points of trajectory in the phase 
space is not easy to be evaluated. In recent years, the study of an increasing 
numbers of natural phenomena that appear to deviate from standard statis- 
tical distributions has kindled interest in alternative formulation of statistical 
mechanics (Kaniadakis 2005). The new formulations of statistical mechan- 
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ics should preserve most of the mathematical and epistemological structure 
of the Boltzmann-Gibbs theory, while reproducing the phenomenology of 
the anomalous systems (Kaniadakis 2005). This K-theory obtained a lot of 
results in the analysis of natural phenomena, but it was not extended to 
irreversible open systems. A modified K-theory has been proposed to use 
it in irreversible open systems, in order to obtain the integral value of the 
Hamiltonian at the end points of trajectory in the phase space, obtaining a 
new statistical expressions for the entropy generation (Lucia 2010). 

The aim of this paper is to find an analytical relation on probability 
to obtain a condition on stability of steady states of open systems which 
represents the condition of quasi-ergodicity for the open systems too. 

To do so in Section 2 the open system is defined, in Section 3 the entropy 
generation is introduced and related to its definition in rational thermody- 
namics, in Section 4 the quasi-ergodicity is discussed in relation to the statis- 
tical properties of entropy generation whose definition related to k— statistics 
is introduced in Section 5 where the required results is obtained. 

2 The dynamical system 

In this section it will be defined the system considered. To do so, we must 
consider the definition of 'system with perfect accessibility', which allows us 
to define both the thermodynamic system and the complex and dynamical 
system. 

Definition 1 (Lucia 2008) - A system with perfect accessibility QpA is a pair 
(i7,n), with VL := {cTj G : (Tj = {pi,q,i),i G [1, A^]}, being N the number 
of particles, e , i G [l,N] their canonical coordinates and Pi G M.^^ , 
i G [1, A^] their conjugate momenta, and U a set whose elements tt are called 
process generators, together with two functions: 



where S is the state transformation induced by tt, whose domain V (tt) and 

range 'TI{'k) are non-em,pty subset ofQ. This assignment of transformation 
to process generators is required to satisfy the following conditions of acces- 
sibility: 

1. IIo" := {Sa : TT G n, a G P (tt)} = , Vo" G Vt: the set Her is called the 
set of the states accessible from a and, consequently, it is the entire 
state space, the phase spase Vt; 





(2) 
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2. if V{7r") n7^(7^') ^ ^ ©(tt'V) = S;,^ {V {tt")) and S^n^,a = 
5^//5^/(7,V(7 e P(7r'V) 

Definition 2 (Lucia 2001) - A process in flpA is a pair (tt, a), with a a state 
and TT a process generator such that a is in V (tt) . The set of all processes of 
QpA is denoted by: 

UoQ = {{7r,a) : TT ell,a eV{7r)} (3) 

If (tt, (t) is a process, then a is called the initial state for (tt, a) and Sa is 
called the final state for (tt, a). 

Definition 3 (Lucia 2001) - A thermodynamic system is a system with per- 
fect accessibility ^pa with two actions W (tt, ct) G R and H (tt, a) e called 
work done and heat gained by the system during the process (tt, cr); respec- 
tively. 

Hypothesis 1 (Lucia 2008) - There exists a statistics [xpa describing the 
asymptotic behaviour of almost all initial data in perfect accessibility phase 
space such that, except for a volume zero set of initial data a , it will be: 

for all continuous functions </? on Q^pa and for every transformation a i->- 

Definition 4 - The triple {flpA,J^,(J^pA), with an algebra or a a— algebra, 
is a measure space, the Kolmogorov probability space F. 

Definition 5 - A dynamical law r is a group of meausure-preserving auto- 
morphisms S : VLpA ^ ^pA of the probability space V . 

Definition 6 (Lucia 2008) - A dynamical system — {^pa, A*pa, t) con- 
sists of a dynamical law r on the probability space V. 

3 Entropy generation 

In thermodynamics, the energy lost for irreversible processes is evaluated by 
the first and second law of thermodynamics for the open systems. So the 
following definition can be introduced: 
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Definition 7 (Lucia 2009) - The entropy generation Sg is defined as: 



Wiost a . Qrf, Ta\ , AH AEk + AEg-W 



1-^ +^-A^+ — ' "^^ (5) 

ref J- a \ J-r / J- a -'-a 

where Wiost is the work lost for irreversibility, Tref the temperature of the 
lower reservoir, Q,. is the heat source, Tr its temperature, Ta is the ambient 
temperature, H is the enthalpy, S is the entropy, Ek is the kinetic energy. 
Eg is the gravitational one and W is the work. 

Following Truesdell (Truesdell 1970), for each continuum thermodynamic 
system (isolated or closed or open), in which it is possible to identify a ther- 
modynamics subsystem with elementary mass dm and elementary volume 
dV — dm/p, with p mass density (Truesdell 1970, Lucia 1995, Lucia 2001, 
Lucia 2008) the thermodynamic description can be developed by referring to 

the generalized coordinates \^i,^i,t> , with the exten- 

sive thermodynamic quantities and af^ their values at the stable states. 

Theorem 1 (Lucia 1995) - The termodynamic Lagrangian can be obtained 
as: 

C = -Tref Sg (6) 

Theorem 2 - The principle of maximum entropy generation (Lu- 
cia 1995, Lucia 2008): The condition of stability for the open system' sta- 
tionary states is that its entropy generation Sg reaches its maximum: 



4 Irreversible ergodicity 

A real-valued random variable vr : QpA — )■ M on F induces a probability 
measurement Ppa '■ ^ [0,1] on the state space (J-']g,R). Considering the 
probability as a property of the generating conditions of a sequence, random- 
ness can be related to predictabihty and retrodictability (Primas 1999). A 
family (t) : t E TZ} is called a stocastic process, which can be represented 
by a family {7 {c (t)) : t G M} of equivalent classes of random variables C, (t) 
on F. The point function {a (t)) is called trajectory of the stocastic pro- 
cess i{t). The description of physical systems in terms of a trajectory of 
a stochastic process corresponds to a point dynamics, while its description 
in terms of equivalent classes of trajectories and an associated probability 
measure corresponds to an ensemble dynamics (Primas 1999). 

Definition 8-^4 stochastic process is said weakly stationary if: 
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1- [e^WL.VteR 

A consequence of the Wiener-Khintchin-Einstein theorem (Primas 1999) 
consists in the existence of a complex- valued function ^ : R — > C, continuous 
at the origin, which is the covariance function of a complex-valued second- 
order, weakly stationary and continuous stochastic process if and only if it 
can be represented as: 

iit)^ r e^^'di{X) (7) 



where, as a consequence of the Bochner-Cramr representation theorem (Pri- 
mas 1999), ^ : R — > R is a real, never decreasing and bounded function, 

called spectral distribution function of the stochastic process. Moreover, the 
Lebesgue's decomposition theorem states that the spectral distribution func- 
tion can be decomposed uniquely as: 

I = c^l<^ + c't + c"'r' (8) 

with > 0, c* > 0, c"'^ > 0, c'^ -h c"* + c"'^ = 1, i'^ step function, continuous 
and singular real function, ^""^ absolutely continuous real function. 

There exists a close relationship between regular and stochastic processes 
and the irreversibility. A system is called irreversible if the lost energy is 
strictly positive. Following Konig and Tobcrgtc (Primas 1999) a linear input- 
output system behaves irreversible if and only if the associated distribution 
function fulfills the Wiener-Krein criterion for the spectral density of a linear 
regular stochastic process: a weakly stationary stochastic process ^ with 
mean value ^ev — and spectral distribution A i— )■ { (A) is linearly regular if 
and only if its spectral distribution function is absolutely continuous and if: 



J — ( 



oo In (A) /d\ 



dX > -oo (9) 



1 + A2 

Following Wiener and Akutowics, we can assume the following: 

Hypothesis 2 (Primas 1999) Every stationary process with absolutely con- 
tinuous spectral function {— regular process) is ergodic. 

The mathematical foundation of ergodic theory is to establish a connec- 
tion between phase average and time average (van Lith 2001). It follows an 
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association between ergodic theory and objective interpretation of probabil- 
ity because the ergodic theory estabhshes a connection between probabihty 
measures and objective features of real word (van Lith 2001). A probabil- 
ity distribution is stationary if it is constant at all fixed points in F, and 
it reflects the fact that the system is in a steady state. Modern ergodicity 
is founded on the concept of measure preserving dynamical systems and 
on the Birkhoff's ergodic theorem (van Lith 2001). Here, considering the 
previous Section, we introduce the stochastic processes. To do so, we state 
that: 

Theorem 3 - If /ipa (^^pa) is finite, then for any integrable function (p : 
^IpA the time average (y?)^ on 7 is defined for all orbits 7 outside 

of a set Af^ of measure jii^M^) = 0. Furthermore {ip)^ is integrable, with 
{(p)^ 07 = {ip)^ wherever it is defined, and with 



CoroUairy 1 - Birkoff's ergodic theorem. A measure preserving trans- 
formation S : VlpA — )■ ^pA is ergodic if and only if, for every integrable 
function 93 : fip^ — )■ M, the time average {ip)^ on 7 is equal to the space aver- 
age ip (a) jipA {da) for all points a outside of some subset Af^ of measure 

Consequently, measurements results as the infinite time averages of phase 
functions because they take a long time compared to the microscopic relax- 
ation time and, for metrically transitive (= ergodic) systems, measurement 
results are almost always equal to microcanonical averages (van Lith 2001). 
Moreover, the initial states lead to different paths in phase space, so the 
averages depend on the initial state. 

Definition 9 - e-steady state. LetV^ be a dynamical system ande = {es} 
be fixed and non zero. An open system is in e— steady state during the time 
interval T if and only if for all S & D,pa there exists ^5 £ K such that for all 
t & T it follows: 



Consequently, the probability may fluctuate within small bounds and, 
consequently, dynamical evolution towards steady states is allowed. Every 
statistical observable is induced by a random variable, while an observable, 
that is a cr-homomorphism, deflnes only an equivalence class of random vari- 
ables which induce this homomorphism. For a statistical description it is not 




(10) 



\{S)r-Cs\<es 



(11) 
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necessary to know the point function a t— t- 5 (a), but it is sufficient to know 
the observable ^. The description of a physical system in terms of an indi- 
vidual function n : a y-> S (a) distinguish between different points a E Qpa 
and corresponds to an individual description of equivalence classes of random 
variables does not distinguishes between different points and corresponds to 
a statistical ensemble description. 

Definition 10 - Let it be n : Qp^ ^ a real-valued Borel function such 
that a 1—^ S {a) is integrable over QpA with respect to fipA, the expectation 
value Tiev of S (a) with respect to fj,pA is: 



As a consequence of the theorem ([2]) and the relation fll5l) . it follows that: 

Theorem 4 - Irreversible ergodicity (Lucia 2008) - In non equilibrium 
transformation the volume of the phase space ^pa contracts indefinitely. 

5 Entropy generation and its statistics 

If the system is ergodic, the time-averages over infinite times coincide with 
the Gibbs phase- averages, but there remains open the problem that noth- 
ing is known concerning the time-averages for large but finite times. It was 
shown that the use of time-averages amounts to introducing a measure in 
phase space, suitably defined by the dynamics of the system. One can imag- 
ine that, if one has to deal with a metastable state ergodic behaviour is 
granted only on a times-scale much larger than the available one, then the 
orbits could exhibit some strange features, such as a non integer fractal di- 
mension on the observed time-scale, which may prevent the use of the Gibbs 
measure. It has been suggested that in some problems where metastable 
states show up, one has to replace the Gibbs measure by the Tsallis one. 
In the present paper we show that, if a system has dynamical time- averages 
compatible with a Tsallis ensemble, then, on a certain time-scale, the orbits 
have a definite non integer fractal dimension. We also show that the diffusion 
process of the orbits in phase space is, in some sense, slower than in the full 
chaotic case corresponding to Gibbs measure. In order to establish a cor- 
respondence between dynamics and Tsallis distribution one has to solve an 
analytical problem by an exact correspondence between dynamics and Tsallis 
distribution, given only for continuous-time dynamical systems and not for 
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mappings. In the case of a mapping, such a correspondence is obtained by 
introducing a suitable Umiting procedure (Carati 2008). 

In 1988 TsaUis postulated the physical relevance of one-parameter gener- 
alization of the entropy (Tsallis 1988), introducing a generalization of stan- 
dard thermodynamics and of Boltzmann-Gibbs statistical mechanics. In 
2002, Kaniadakis proved some properties for the generahzed logarithm, so 
that we can introduce the following: 

Definition 11 (Kaniadakis 2003, Pistone 2009) The generalized exponential 
exp^^} is defined as follows: 

exp,,,(r) = exp ( f = [ (- + v^TT^)^^^ ^ 

(13) 

Definition 12 (Kaniadakis 2003, Pistone 2009) The generalized logarithm 
Inj^} is defined as follows: 

I-W(^) - (14) 

Definition 13 (Lucia 2010) The statistical expression, for the irreversible- 
entropy variation, results: 

5g = - J^p^ln{„}(p^) (15) 

7 

with 

p^ = a exp^^y ^ - ^\rp^ ^ (16) 

p^ is the probability of the path 7, and a — [ exp ( — 

i/y —^^'^^^^^^dV + ■^)] ^ , with H = u — J2il^iPi non- equilibrium gen- 
eralisation of the grand- canonical Hamiltonian, H^(q-^ and H^(^^^ the values of 
H at the end points of trajectory 7 m r = and t, and the time-averaged 
rate entropy production ofj, and A = 2kB are two arbitrary, real and positive 
parameters, /i — ^1/^2, with (3i and ^2 Lagrange multipliers in the entropic 
functional 

^=^«-/3i||x^p^-ij -p^l^Y.P'^^^-^^ (17) 
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which must be stationary for variations of jSi and (32, and where is the 
energy of the j-th path and U is the total energy, ln{K}(x) is an arbitrary 
function, a generalization of the logarithm function, whose properties have 
been studied by Csiszar (Csiszar 1967) and when ln^i^j{x) = \n{x) the gener- 
alized entropy reduces to the Boltzmann-Gibbs- Shannon entropy (Wada and 
Scarfone 2009). 

Lemma 1 Let us consider the In^p^. Its first order derivative is: 



^7 



Proof. 

From the Definition [T^ it follows that: 



K(P7) 




(19) 



2 k p^ 2k 



K 



ln.(p^)+j9;('^+i) 

P7 



□ 



As a consequence of the theorem of maximum entropy generation the 
following inequality can be obtained the following condition: 

Theorem 5 The condition of stability for the open system ' stationary states 
is the following inequality for the probability of the paths: 

-E.lnw^7^^ (20) 



Proof. 
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From the Definition [T3] it is possible to write: 



(21) 



7 7 V / 

7 7 7 

7 7 
and, considering the Theorem [21 it follows: 

S'^>0^-Y,i^ + ^) ln.(p,) - E^;' > (22) 

7 7 

that allows us to obtain the following condition of stability: 



> T— (23) 



□ 

6 Comments and conclusions 

As Berkovitz, Frigg and Kronz pointed out (Berkovitz et al. 2006): 'The 
relevance of ergodic theory for the effective modelling of physical systems 
has come into serious question due to KAM theory. Kolmogorov (1954) for- 
mulated the theorem, which concerns Hamiltonian systems and perturbation 
theory (...) The theorem basically says that if a small perturbation is added 
to an integrable system with two physical degrees of freedom, then the tori 
with a sufficiently irrational winding number survive the perturbation. (...) 
Markus and Meyer (1974) developed some important consequences of the 
KAM theorem (...) The Markus- Meyer theorem is often informally char- 
acterized as the claim that generic Hamiltonian dynamical systems are not 
ergodic. This is a striking claim and its effect appears to undermine the 
relevancy of ergodic theory immediately. But this effect is mitigated as soon 
as the theorem is stated in full: the set of Hamiltonian flows that are ergodic 
on the energy hypersurfaces associated with each element of a dense set of 
energy values is of flrst category in the set of inflnitely differentiable Hamil- 
tonian flows. In short, the class of ergodic Hamiltonians is of flrst category in 
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the set of generic Hamiltonians. A set P is of first category in set Q,if P can 
be represented as a countable union of nowhere dense sets in P. All other 
sets are of second category. Loosely speaking, a set of first category is the 
topological counterpart to a set of measure zero in measure theory (...) The 
contrast class is second category, which is (. . . ) is the topological analogue 
of a set of non-zero measure (...) The key phrase that has the mitigating 
effect is infinitely dijferentiable. This restriction is substantial. It actually 
rules out a large class of Hamiltonian systems that are of physical interest 
particularly for classical statistical mechanics (...) it is worth contrasting 
the Markus-Meyer theorem with a theorem (the first) of Oxtoby and Ulam 
(1941): in the set of measure-preserving generahzed dynamical flows, the 
set of dynamical flows that are not ergodic is of flrst category. The class of 
Hamiltonian flows of physical interest is substantially broader than the class 
considered in the Markus-Meyer theorem, but substantially narrower than 
that considered in the Oxtoby-Ulam theorem. (...) Boltzmann (1871) stud- 
ied hard ball systems in developing a mathematical foundation for statistical 
mechanics. He conjectured that such systems are ergodic when the number 

of hard balls is sufficiently large. Whether this is the case is still an open 
question. However, it turns out that hard ball systems are ergodic for some 
small A^. Sinai (1970) provided the first rigorous demonstration that a hard 
ball system is ergodic. (...) the Boltzmann-Sinai ergodic hypothesis (Sinai, 
1963) (...) says that a system of N hard balls on or is ergodic for 
any N > 2, and it was inspired by Krylov's observation, made in 1942 that 
hard ball systems exhibit an instability that is similar to geodesic flows on 
hyperbolic surfaces, known at the time to be ergodic.'. 

Ergodicity is the mechanism that allows us to pass from a mechanical 
description of the system to a statistical one. A fundamental question is 
what happens if the system is not ergodic. If the system is not ergodic, we 
can still have a relatively simple statistical description of it if the trajectory 
is a fractal curve embedded in the smooth energy surface of fl cells. If 
the reduced box-counting dimension of such a curve is < d < 1 then 
the available phase space is a fractal object embedded in the Boltzmann 
6A^— dimensional smooth phase space: when d = 1 we regain completely the 
ergodicity (Garcia-Morales and Pellicer 2006). Starting from the analysis of 
the entropy and information on fractal space Naschie (Naschie 1992) showed 
a connection between the information and the topological dimension of some 
fractal sets. Moreover, he showed that the dynamics of this kind of phase 
space is quasi-ergodic for a dimension A^ > 4 of the number of particles in 
the fractal phase space of topological dimension df (Wang 2004). 

The unifying principle of the maximum entropy generation allows to build 
a generalized thermodynamic formalism analogously to that of Hamiltonian 
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mechanics. This is clearly seen in the definition of the entropy (jH]) and f[T^ . 
Many experiments in various fields of nuclear and condensed matter physics 
suggest the inadequacy of the Boltzmann-Gibbs statistics and the need for 
the introduction of new statistics (Kaniadakis and Lissia 2003). In particular 
a large class of phenomena show power-law distributions with asymptotic 
long tails. Typically, these phenomena arise in presence of long-range forces, 
long-range memory, in chaotic systems, in opens systems or when dynamics 
evolves in a non-euclidean multi-fractal space-time. In the last decade many 
authors pursued new statistical mechanics theories that mimic the structure 
of the Boltzmann-Gibbs theory, while capturing the emerging experimental 
anomalous behaviors (Kaniadakis and Lissia 2003). In the last ten years, 
an intense activity has been modifying and improving our understanding of 
contemporary Statistical Mechanics and Thermodynamics. 

The aim of this paper was to find an analytical relation on probability 
to obtain a condition on stability of steady states of open systems. This 
is the relation ( l20l) . It represents the condition of quasi-ergodicity for the 
open systems and it is a consequence of the principle of maximum entropy 
generation. 
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